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Abstract. In a statistical study of the decay of individual sunspots based on DPR data 
we find that the mean instantaneous area decay rate is related to the spot radius r and the 
maximal radius ro as D = Cd r/ro, Cd = 32.0±0.26 MSH/day. This implies that sunspots 
on the mean follow a parabolic decay law; the traditional linear decay law is excluded by 
the data. The validity of the Gnevyshev-Waldmeier relationship between the maximal 
area Aq and lifetime T of a spot group, Ao/T ~ 10 MSH/day, is also demonstrated for 
individual sunspots. No evidence is found for a supposed supergranular “quantization” of 
sunspot areas. Our results strongly support the recent turbulent erosion model of sunspot 
decay while all other models are excluded. 

Abbreviations: CM - central meridian; DPR Debrecen Photoheliographic Results; 
GPR - Greenwich Photoheliographic Results; MSH - millionth solar hemisphere; MSHER 
- MSH equivalent radius 


1. Introduction 

It has been known since the era of the first telescopic observations that 
sunspots disappear from the face of the Sun on a timescale of days/weeks 
by a gradual shrinking (sometimes punctuated by fragmentation events). 
During this process of decay the boundary of the spot remains well defined 
and its temperature deficit does not change too much (nor does its mag¬ 
netic field strength vary by more than about 20%, cf. Collados et al, 1994; 
Steinegger et al., 1996). The extremely large individual variations among 
spots however did not make it possible to find quantitative regularities in 
the decay until the availability of sufficiently large and precise observation¬ 
al databases (GPR, Pulkovo, Mt. Wilson). These studies began in earnest 
towards the middle of this century, and by about 1970, two simple statistical 
relationships had become rather widely accepted among solar physicists: 
The rule of the proportionality of the maximal area Aq of a sunspot 
group to its lifetime T (first plotted by Gnevyshev, 1938, and formulated 
by Waldmeier, 1955): 

Aq = DqyjT Pgw ~ 10 MSH/day. 


(1) 
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The alleged linear time-dependence of the area during the phase of 
decay: 

A = Ao — D(t — to). (2) 

This formula was originally proposed by Bumba (1963) for recurrent 
sunspot groups, while for the majority of non-recurrent groups a concave 
decay curve was found to be a better fit. However, for a variety of 
reasons including bad sampling intervals and large scatter, significant 
departures of the decay curve from linear were difficult to find and even 
more difficult to describe quantitatively. Thus, the linear decay law had 
become widely accepted as the simplest approximation to the decay 
curves (Ringnes, 1964a, 1964b; Robinson and Boice, 1982). 

More recently, doubt has been cast on the validity of the linear decay law. 
Moreno-Insertis and Vazquez (1988) and Martinez Pillet, Moreno-Insertis, 
and Vazquez (1993) pointed out that the statistical evidence for a parabolic 
decay law 


A = Aq — 2 \J TrAow(t — to) + Tiw 2 {t — to) 2 (3) 

is somewhat stronger than for a linear decay. Such a decay law would follow 
from a slow inwards motion of the spot boundary with a constant velocity 
w. However, a definite conclusion about the significance of that result could 
not be drawn. The correct form of the sunspot decay law has thus remained 
undetermined. 

Note also that the Gnevyshev-Waldmeier relation (1) was derived for 
sunspot groups, and its validity for individual spots has not been demon¬ 
strated. 

Another open question is that of the existence of preferred length/area 
scales corresponding to higher stability of the spots (i.e. to lower decay 
rates). Such a “quantized” character of the distribution of sunspot areas 
(possibly caused by interaction with supergranulation) has been proposed 
by several authors (Ikhsanov, 1967; Dmitrieva, Kopecky, and Kuklin, 1968; 
Bumba, Ranzinger, and Suda, 1973). The statistical significance of those 
findings however remains dubious. 

In this paper we present the results of a statistical investigation aimed 
at resolving these fundamental open questions related to sunspot decay. 
(Preliminary results of this study have been published elsewhere: Petrovay 
and van Driel-Gesztelyi, 1997.) Our investigation yields definitive answers 
to the issues mentioned above. This has been made possible by our use of 
the Debrecen Photoheliographic Results: this recently published database 
is the only data set of statistically significant size to contain parameters for 
individual sunspots with day-to-day identification, measured and reduced 
by very high standards of precision. 
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The available theories and models of sunspot decay served as our guide 
during the analysis by concentrating our attention to those correlations for 
which different models had different predictions. We are thus also able to 
exclude some of those models, leaving a single one compatible with the 
observational constraints. The models and their predictions are outlined in 
Section 2. In Section 3 we discuss our selection and handling of data in 
detail, while Sections 4 and 5 present the results. Section 6 concludes the 
paper. 


2. Models 

The first quantitative model for sunspot decay was proposed by Gokhale 
and Zwaan (1972). To explain the observed sharp boundary of spots they 
assumed that the magnetic flux tube is surrounded by a current sheet, i.e. 
that moving outwards from the spot center, the magnetic flux density drops 
to zero in a very narrow radius range. They further assumed that turbulence 
is fully inhibited by the strong magnetic field in and inside the current sheet, 
so that the decay of the spot is only due to pure Ohmic dissipation in the 
current sheet. They pointed out that a linear decay law implies that the 
thickness of the current sheet must be a constant fraction of the radius, 
i.e. the spot remains self-similar during the decay. In order to produce an 
ever thinner current sheet, however, they needed to postulate ad hoc the 
existence of a radial inflow towards the tube, for which no observational 
indication existed. For this reason, their model never became so widely 
accepted as the turbulent diffusion model. 

The latter model, proposed independently by Meyer et al. (1974) and by 
Krause and Rudiger (1975), assumes that the decay of the spot is caused by 
turbulent diffusion throughout the whole cross-section of the tube, with a 
magnetically reduced but otherwise constant diffusivity. It was found that 
defining the boundary of the spot at some arbitrary fixed value of the flux 
density, such a diffusion naturally results in a linear decay law, and the 
Gnevyshev-Waldmeier relation is also returned. The diffusivity was essen¬ 
tially treated as an adjustable parameter. The tubes in this model had no 
sharp boundary, being gradually “washed away” by diffusion, and the cen¬ 
tral flux density showed a gradual decrease instead of remaining more or 
less constant until the full decay of the spot. The model is also problematic 
from a theoretical point of view, as a constant magnetically reduced dif¬ 
fusivity is hardly compatible with a magnetic field that must obviously be 
lower outside the tube. Despite these obvious shortcomings, its prediction 
of the two (alleged) linear laws of sunspot decay made this model the most 
accepted theoretical explanation of the decay for two decades. 


decaylpp.tex; 1/02/2008; 4:42; no v.; p.3 



4 


PETROVAY AND VAN DRIEL-GESZTELYI 


The studies of Moreno-Insertis and Vazquez (1988) and Martinez Pillet, 
Moreno-Insertis, and Vazquez (1993) questioned the validity of the linear 
decay law, reviving the old suggestion (Simon and Leighton, 1964) that the 
decay of a spot may proceed by the eroding action of external turbulence 
that “gnaws” off bits and pieces of the flux tube. It was time to recon¬ 
sider the quantitative decay models. The recent turbulent erosion model 
(Petrovay and Moreno-Insertis, 1997) is formally a generalization of the 
turbulent diffusion model allowing a variable diffusivity explicitly depend¬ 
ing on the magnetic field. The solutions of the nonlinear diffusion equation 
in cylindrical geometry 




f v{B) 


dB 

dr 


( 4 ) 


were investigated by numerical and analytical methods. Here, B is the mag¬ 
netic flux density, r is the radial coordinate, and v is the turbulent magnetic 
diffusivity. For the dependence of the diffusivity on the magnetic field they 
considered a simple function satisfying the basic physical requirements: 

tT>\ u o 

" (B) = 1 + \B/BJp 

Here a u is a parameter quantifying the steepness of the diffusivity cutoff 
near B e , the latter being the field strength where the diffusivity is reduced 
by 50%. Physically, one expects B e ~ B eq , the kinetic energy density of 
photospheric turbulence. The solutions of Equation ([|) were found to be 
qualitatively similar to the case with constant diffusivity as long as the 
inhibition of turbulence was weak (i.e. for low values of a v and/or Bo/B e , 
Bq being the central field strength of the initial field profile). For strong inhi¬ 
bition of turbulence however a new class of solutions sets in: a current sheet 
is spontaneously formed around the tube, moving inwards with a constant 
speed w, while the field strength inside the current sheet remains practically 
unchanged during the decay until the arrival of the current sheet. This leads 
to a parabolic decay of the cross section of the tube. 

An important and attractive feature of the erosion model is that the 
velocity w of the current sheet becomes asymptotically independent of a u , 
i.e. for a v —> oo it goes to the finite limit 


( 5 ) 


VJ 


0-1/3 -^e *4) 

H 0 - B e r 0 ' 


( 6 ) 


where ro is the maximal radius of the spot. Thus, the lifetime of the spot 
remains finite even for infinitely effective inhibition of the turbulence inside 
the tube. 

Using Equation (|) and assuming that the unperturbed value of the 
diffusivity outside the tube is uq = 1000km 2 /s, i.e. the granular value, 
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Self-similar 

sunspot 

Turbulent 

diffusion 

Decay law 

D = const. 

D = const 

Lifetime 

T oc Ao 

T oc To 

B 0 (t) 

constant 

gradual decrease 

Spot boundary 

sharp 

diffuse 

Reference 

Gokhale & 

Meyer et al. (1974) 


Zwaan (1972) 

Krause & Rudiger (1975) 



Turbulent 

erosion 

Universal 

parabolic 

Decay law 

D oc r/ro 

D oc r 

Lifetime 

T oc Ao 

T oc \A4o 

B 0 (t) 

constant 

constant 

Spot boundary 

sharp 

sharp 

Reference 

Petrovay & 

Martinez Pillet 


Moreno-Insertis (1997) 

et al. (1993) 


£>o = 3000 G and B e = 400 G, one finds that the Gnevyshev-Waldmeier 
rule is returned: 

ro/w ~ 2 1 / 3 ——— ——■ = Aq jD D ~ 10. (7) 

B e Vo 

For the instantaneous decay rate D = A in turn we have the prediction 

D ~ (V^ttuq Be ) — = C D - C D ~ 22. (8) 

V Bo - B e J r 0 r 0 

The predictions of different published decay models are summarized in 
Table I. The “universal parabolic” model in the table is identical to the 
turbulent erosion model except that the current sheet velocity w is assumed 
to be a universal constant instead of depending on the maximal spot radius 
as in Equation (|6|). This “model” has no theoretical foundation, and it was 
only tentatively suggested by Martinez Pillet, Moreno-Insertis, and Vazquez 
(1993) to forge a link between the lognormal distributions of decay rates and 
spot areas. 
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3. Data 


3.1. Selection 

DPR data for the years 1977 and 78 (Dezso, Gerlei, and Kovacs, 1987; 
Dezso, Gerlei, and Kovacs, 1997) were used, selecting only spots with max¬ 
imal areas exceeding 10 MSH and measurements made within ±72° of the 
central meridian. In the following, this part of the solar surface will be 
referred to as “the visible hemisphere”. For the present study only total 
(umbra+penumbra) areas are used. Altogether, our data set consisted of 
3990 area measurements for 476 different spots. 

The nominal error of measurements in the DPR is less than 1 MSH, 
implying a “worst case” visibility function (Kopecky, Kuklin, and Starkova, 
1985) sec 2 A where A is the angular distance from the center of the disk. 
Variable exposure times and seeing conditions, aggravated by the geometric 
foreshortening correction factors for higher A’s however lead to a larger ran¬ 
dom error in the area measurements. The Introduction of the DPR (Dezso, 
Gerlei, and Kovacs, 1987) quotes error values of typically 5 %, reaching 20 % 
for the smallest spots in the sample. In spite of the great precision and care 
the DPR was compiled with, there is also a small chance for errors in the spot 
identification, especially in complicated, quickly changing sunspot groups. 
These errors are small compared to the random physical scatter present in 
the data, and being random, they cannot contribute to the physical corre¬ 
lations we find in Sections 4 and 5. 

A systematic error is also involved owing to the physical foreshortening 
(Archenhold, 1940). While this error can be quite significant for small spots 
close to the ±72° CM distance limits, it will not exert a great influence on a 
sample as a whole unless the spots of that sample have a strong tendency to 
be situated close to the limits. Such a tendency is known to be present in at 
least one case. The first observation of a spot on the “visible hemisphere” is 
considered as its birth if the spot was not seen at the previous observing time 
(typically, the previous day). The death of a spot is defined in an analogous 
manner. While the distribution of birth and death sites has a well-known 
East-West asymmetry, peaking near the limbs, a straightforward estimate 
following the method of Kopecky, Kuklin, and Starkova (1985) shows that 
with the visibility function quoted above the effect is small, and birth and 
death sites in our sample are distributed nearly uniformly, i.e. the systematic 
error is small. 

One should also be aware that some steps in the compilation of the data 
set lead to selection effects. An obvious effect is that the DPR database only 
includes spots for which at least 2 observations exist. Taking into account 
that observations are typically separated by 1-day (and occasionally 2 to 4- 
day) intervals, this implies that our sample is not complete below lifetimes 
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of about 5 days. In what follows, this effect will be referred to as the short 
lifetime selection effect. 

Another effect is the visibility reduction of lifetime (Kopecky, Kuklin, 
and Starkova, 1985) which will cause a systematic error downwards in deter¬ 
mined spot lifetimes. The error is maximal for lifetimes of 8-10 days, but 
with the visibility function at hand it altogether remains small. 

3.2. Definitions 


Let us assume that a given spot was observed at times ii,... ,tjv, the cor¬ 
responding area values being A\,.... A^. If t\ and tjsr correspond to the 
actual birth and death of the spot, as defined in Section 3.1 above, then the 
lifetime of the spot is given by 

T d = t N -t i T l d (9) 


For each measured area value the equivalent radius r* of the circle having 
the same area is computed. In this paper, equivalent radii will be expressed 
in MSHER units: a spot of radius 1 MSHER has an area of 1 MSH, so 
that in these units A, L = r'f (1 MSH= 3044km 2 , 1 MSHER= 984km.) Time 
will be measured in days. The spot is regarded to be in a decay phase at 
time ti if i > > rt + 1 holds (or if rjy-i > r at, in the case of the last 
observation). The “instantaneous” decay rate of the spot is then given by 
Di = 27 xriWi where w is the velocity of the spot boundary, calculated from 


w, = + 2 K lll -d° > )fe -«■:-!) 


(b+1 ti— 1)“ 


w {j) _ D-l+j fi+i 
ti+j b— i +j 


( 10 ) 

( 11 ) 


In the case of the last observation, w is calculated for a time tjv' = (bv-i + 
tjsr)/ 2 from 


W N ' 


r N—i - r N 
tN ~tN -1 


( 12 ) 


We note that the above formulae imply the assumption w = const, between 
observations, i.e. a parabolic decay law. In order to ensure that our results 
are not distorted by such a tacit assumption, we repeated all statistical 
studies by using an alternative definition for D, based on a linear decay 
law. The results of Sections 4 and 5 were not significantly influenced by this 
choice. As those results exclude the linear decay law, here we use the above 
formulae in all calculations. 

Similar care must be exercised with the definition of the maximal spot 
area Aq, corresponding to an equivalent radius ro, at time to- Many complex 
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spots show several local maxima (and consequently several decay phases) in 
their A(t) curves, so the definition of a single maximum is not unique. For 
the results presented in this paper the maximum was defined as the absolute 
maximum preceding the given observation of the spot: 

= Max {Aj}j=i,...,i Aq = Aq,n- (13) 

Again, using alternative definitions for the maximum, the results were 
not significantly influenced (although using the preceding local maximum 
instead of the absolute maximum made most correlations somewhat worse). 

Another problem related to maximum determination is that owing to the 
rotation of the Sun, in many cases we cannot follow the complete evolution 
of the spot. In these cases we cannot be sure if the absolute maximum we 
determined is really the absolute maximal area of the spot during its prior 
evolution. In order to deal with this problem, the sample was divided into 
subsamples of different “reliability” (e.g. spots born and died on the visible 
hemisphere etc.). While using a more carefully selected sample does reduce 
the scatter in statistical relationships, it may also reduce the significance of 
the result owing to the lower number of data. Again, most of the studies 
presented in Sections 4 and 5 were performed on different subsamples. The 
results presented here use subsamples chosen to minimize the scatter while 
being still large enough to be statistically significant. 

3.3. Bias corrections 

The explicit selection effects mentioned in Section 3.1 above originate from 
conscious decisions during the compilation of the data set. Our sample is 
however also affected by implicit selection effects or biases. A bias here 
means that the distribution of some variable y in our data set will differ 
from its actual distribution among sunspots because the expected number 
of observations for one spot depends on y. Let P{y; x) denote the observed 
distribution of y in some subset of our data homogeneous in the variable(s) 
symbolically denoted here by x. If the real distribution (i.e. the “spot dis¬ 
tribution” as opposed to the “data distribution”) is denoted by P(y;x) and 
the expected number of data for one spot by p(y; x) then clearly 

P(y, x) = F(y ; x)P(y; x) (14) 

where F oc 1/p, with some arbitrary normalization. It is convenient to 
choose the normalization so that in our discrete sample of N data, F.\ = 

N. Formulae for p viz. for the correction factor F for cases of practical 
importance will be given below. 
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Often we are only interested in the mean and the standard deviation of 
P. These and their r.rn.s. errors are given by 


N 


y = ( yP) = -^ E yiFt 


i=l 


a 2 = ((y-y) 2 P) = 


1 


N 


O- 


cr 


= " 2 E 


("-Dfe 

W /a-x2 a 2 N 

^ \dy. 

8a\ 


J2(yi~y) 2 Fi 


2 2 \ ~~ i ^y _ \ 2 _ a “ \ ' (i 7 i dFj \ 

- a E ( “ Z\r 2 2 ^ (- ' 


i =1 


dyiJ 


- = 


^ \dyi) 

N r 

E 


i 


4(iV-l) 2 ^ 


OFF 2 


2{yi -y)Fi + (yi - y) 2 -^- 1 

oy 


(15) 

(16) 
(17) 


(18) 


If one wishes to make a least-squares fit y = yo( x ) to the data, the 
simplest way to do this for a biased sample is to bin the data in x using 
bins small enough that yo may be securely assumed to be nearly constant 
within each bin. Then y is calculated for each bin from Equation (15), and 
the fit is performed on these bin averages, weighting each value by the 
inverse of its cr|, as given by Equation m- 


3.3.1. y = f(T): Gnevyshev-Ringnes correction 

The lifetime of a spot can obviously only be determined if both its birth 
and death occur in the visible hemisphere, so only such spots will enter our 
data set. The probability p of this being so (i.e. the expected number of data 
for one spot) however depends on T according to a formula first given by 
Gnevyshev (1938) and later corrected by Ringnes (1964c). For our purposes 
this can be written as 


0 


if \a\ > 28 0 


2#o — | a 

360° 


if |a| < 20 q 


(19) 


a == [180° - (SIT + 180°) mod 360°] 


( 20 ) 


Here, $o = 72° is the CM distance limit and H = 360°/P S y n , -P sy n being the 
synodic rotation period of the given spot. For simplicity, here we use the 
value P syn = 28 d for all spots. (In 1977/78 sunspots were situated at high 
latitudes.) 
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The correction factor F is the normalized inverse of Equation (|l9|) , inde¬ 
pendently of the form of the function f(T). Its derivative F'(y), entering 
Equations ( |17P and (18), will obviously depend on that form, though. In 
Section 5 below we will deal with the distribution of y = log T: in that case 
F'(logT) = lnlOTF'(T). 


3.3.2. y = f[D;g(A,Ao)\: correction for the decay law 
Let us first restrict our interest to the case g(A, Aq) = A. In this case our 
study of the distribution of y is limited to spots in the area interval [A, A + 
dA], A spot with a decay rate D will spend a time dA/D in this interval, 
so if our observations are made with l d intervals, the expected number of 
observations for one spot in that interval (in our units of measurement) is 
clearly p = dA/D. In the more general case 


dA dg 

P= dg^' 


( 21 ) 


It is apparent that p (and F) here depend on the form of the g(A,Ao) 
function. This is in contrast to the Gnevyshev-Ringnes correction where 
the observation probability of a spot of lifetime T was independent of any 
other spot parameter than T. For the two cases investigated in Section 4 
below the correction factor takes the form 


’ if g = (AMo) 1/2 

A 0 g 

F 0 — if g = A 1 / 2 

g 


( 22 ) 


Fq being a normalization factor. F is again independent of the form of f(D) 
while its derivative for y = log D is F' = In 10E. 


3.3.3. Limitations of the bias correction 

Assume that all our measured yi data fall within an interval J. If p> 1/N 
in the neighbouring intervals of similar size then we can be reasonably sure 
that the lack of data points outside the above range is not an artefact caused 
by the bias. (This is the case with our studies in Section 4.) If however the 
range of strong bias p ^ 1/N overlaps J or is immediately adjacent to it, we 
should suspect that the lack of data in the strong bias range is an artefact. 
In this case, no proper correction is possible. The values derived for the 
mean and standard deviation, however, may still be acceptable [albeit with 
the large error bar given by Equation (|l^)] if data points are present both 
above and below the strong bias range, thereby offering the possibility of 
an interpolation. This is the case we will encounter in Section 5. 

Beside the implicit selection effects dealt with above, there is a variety of 
subtle biases affecting certain distributions, the corrections of which is not 
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always straightforward or worth the effort. We will meet two such examples 
in the following section. Unfortunately, neither these more involved selection 
effects nor others discussed in Sections 3.1 and 3.3 were properly taken into 
account in much of the earlier work on the subject. 


4. Decay law 


4.1. D VS. r/r 0 

Guided by the prediction of the turbulent erosion model, Equation (||), 
in Figure |] we plot the instantaneous decay rates D against the relative 
radius r/r 0 [in the sense of rj/r'o^, cf. Equation (j^)]. In this plot we show 
measurements for those spots which died on the visible hemisphere (888 
data points, 170 of these for recurrent spots). The results for the complete 
sample are however very similar. As sunspot decay rates are lognormally 
distributed (Martinez Pillet, Moreno-Insertis, and Vazquez, 1993; see also 
the following paper in this series), for the application of Gaussian-based 
statistical methods (least squares etc.) D should be plotted on a logarithmic 
scale. 

It is apparent that the formula (|j) is indeed a reasonable representation 
of the data. Binning the data in r/ro and correcting the bias described in 
Section 3.3.2 we find that a constant D (i.e. a linear decay law) can be 
excluded at a very high confidence level (Q = T[(IV — l)/2,x 2 /2] < 10” 5 ). 
In contrast, a linear fit of the form D = Cpr/v o (i.e. a parabolic decay 
law, the best fit being Cd = 32.0 ± 0.26) yields an excellent agreement 
(Q = 0.68). The difference in the coefficient compared to Equation @ can 
be easily accommodated taking into account the uncertainty of the value 
of ZA), the strong height-dependence of B e , and the fact that our formula 
(Hj) was based on the infinite ot v limit (]6|), so the actual value of w may be 
somewhat higher. 

The question arises, how is it possible that this obvious nonlinearity of 
the decay had not been noticed earlier? While our use of the DPR with 
its high precision individual sunspot area data may be a factor in a more 
conclusive result, the study of Moreno-Insertis and Vazquez (1988) shows 
that with proper methods of analysis a nonlinearity can also be detected 
in the GPR data. We should therefore suspect that the choice of the tech¬ 
nique of analysis is a crucial factor in this respect. Sunspot decay curves 
are usually charactarized by bad sampling intervals: for short-lived spots 
one often only has 2-3 measurents in the decay phase, while for recurrent 
spots a large number of measurements are grouped together in the vicinity 
of 2-3 points of the decay curve. The true endpoint of the decay curve is 
thus ill determined and the method of producing “typical” decay curves by 
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All spots Recurrent spots 





Figure 1. Decay rates plotted against relative radius for the original (top) and binned 
(bottom) data. Data for spots which died on the visible hemisphere are shown. Solid line: 
Equation (8). Dashed line: linear least-squares fit. (Here and in all other figures, 2a error 
bars for the mean are shown.) 


time-shifting individual curves to a common endpoint (e.g. Ringnes, 1964a, 
1964b) implies rather large errors. These errors are further aggravated by 
the fact that the lifetime-maximal area correlation (cf. Section 5) involves a 
large scatter, thus spots with greatly different tq values will enter the sam¬ 
ple, leading to different decay rates at the same time. The resulting large 
scatter in D will tend to smear out any nonlinearity. This is borne out in 
Figure [| where the parabolic decay curve is directly seen in the subsample 
with Aq in a narrow range (bottom left panel), while the picture is much 
more confuse if no such selection is performed (top left panel). (Note that 
some non-linearity is still visible even in this latter panel. This may be due 
to our time-shifting the curves to a common maximum instead of a com¬ 
mon endpoint, as usual.) The selection of an Alp-limited subsample however 
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Non —recurrent spots 



Recurrent spots 



(t-to) 



Figure 2. Relative area vs. time plots. Top: all spots; bottom: spots with maximal areas in 
a restricted range (90-110 for non-recurrents, 250-350 for recurrents). Solid line: Equation 
(3) with Aq = 90 and w from Eqs. (6) and (8) 


greatly reduces the sample size making it difficult to establish the nonlin¬ 
earity at a high confidence level for all but the most frequent Aq values. 
Thus, the D-r/ro plot (Figure [l|) is much more useful for the study of the 
decay law. 

The study of Bumba (1963) resulted in a nonlinear decay law for the 
majority of non-recurrent spots while a linear decay law was found to 
describe the recurrent spots. To check the reality of such a distinction, in the 
right-hand panels of Figure [I] we present the D-r/r$ plot for recurrent spots 
only. Essentially the same results are found as for the complete sample. This 
is so despite the fact that this plot is distorted by an uncorrected bias: the 
apparent local minimum near r/ro — 0.7 is caused by the fact that between 
the first and second disk passages, higher-H spots shrink more —so the 
corresponding gap in r/ro will be wider in the upper part of the diagram. 
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Bumba’s finding could be due to the fact that the Aq values of recurrent 
spots are scattered within very wide limits. This, together with the bad 
sampling intervals, leads to a very confuse area-time plot (righ hand panels 
in Figure |2|), for which a linear fit is as good as anything else. 



Figure 3. Histogram of maximal radii (a) and r vs. r/ro plot (b) for all spots in our data 
set. 

4.2. D VS. r 

Having excluded the linear decay law and shown that D oc r/ro is compatible 
with the data, we still need to check the validity of the prediction of the 
fourth model in Table I, i.e. the “universal parabolic” decay law D oc r. This 
is necessary because in our data there is a good, nearly linear correlation 
between r and r/ro, caused by the rather sharp peak in the distribution of ro 
values, as illustrated in Figure []. (Note that this histogram is contaminated 
by several selection effects, but this is irrelevant for our present purpose.) 
The suspicion arises that the correlation between D and r/ro found above 
might be just a reflection of a more basic relationship between D and r. 

The distribution of the data points from the full sample (1234 points) 
on the r—D plane is shown in Figure |3]. It turns out that a fit of the form 
D oc r, corresponding to the universal parabolic model, is incompatible with 
the data (Q < 10 -5 ), as is a constant decay rate. 

The actual distribution of the data finds its explanation in Figure [|d. 
For r ^ 7, there is a nearly one-to-one correspondence between r and ro, 
i.e. most of these measurements refer to spots close to their maximal radii. 
This is a consequence of the rather fast decrease of the ro distribution curve, 
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Figure 4- D-r plot for the original (a) and binned (b) data, (c): low r part with higher 
resolution; (d): r vs. ro plot. All decay rate measurements are shown. Solid line: linear 
least-squares fit. 


Figure §b in this regime. According to Equation (||), the decay rate is 
then constant: indeed, for r ^ 7 the D-r plot is compatible with D = 
constant. The sharp peak in the histogram of ro near 10 also implies that 
for lower values of r most of the measurements are for spots with ro ~ 10 = 
const. Thus, as it follows from Equation (|8|), in the range r = 3-7 the D-r 
distribution is reasonably well fit by the linear relation D oc r. 

In the range r ^ 3 the distribution again seems to depart from the linear 
law, this time being steeper. This is the consequence of a subtle bias: the 1- 
day sampling interval is relatively “coarser” for small (and thus short-lived, 
cf. Section 5) spots which are thus less likely to be “caught” in the latest 
stages of their decay. Consequently, the mean value of ro begins to increase 
again with decreasing r in this regime, Figure |||d. 
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The fact that the D—r relation for r > 7 is compatible with a constant 
decay rate implies that we have no significant evidence for the existence 
of a supergranular area “quantum” where spots are supposed to be more 
stable, decaying slower. An independent test of this was also performed by 
taking the power spectrum of the histogram of r (not shown here): again, 
no significant peaks were found. The area histogram published by Howard 
(1996) leads to similar conclusions. While it cannot be excluded that using 
a much larger database, such an effect may be discovered in the future 
(note the slight local minimum in the D-r plot near r ~ 14), none of the 
papers that claimed to have found such a quantization was in fact based on 
a database significantly more extensive than ours in the given area range 
(<; 80MSH). Some of these works used data for “regular” spots only, the 
selection criteria however were rather subjective. Facing this state of the 
matter, we are forced to conclude that no compelling evidence exists for the 
supposed phenomenon of quantized sunspot areas. 



1.6 1.8 2.0 2.2 2.4 2.6 
log A„ 


2.0 



0.0 . 

1.6 1.8 2.0 2.2 2.4 2.6 
log A„ 


Figure 5. Area-lifetime relation for all spots with determined lifetimes in the 1977/78 
DPR data. +: non-recurrent spots; O: recurrent spots. Solid line: Gnevyshev-Waldmeier 
relation Ao = 10 T; dashed line: linear best fit. 


5. Area-lifetime relation 

There are 128 spots in our sample which were born on the visible hemi¬ 
sphere and also died there, so that their lifetime could be determined. These 
lifetimes are plotted against the maximal radii in Figure ||. The empty hor¬ 
izontal band corresponds to spots living for about half a solar rotation: 
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obviously, no lifetime could be determined for any of these. After binning 
and applying a Gnevyshev-Ringnes correction (Section 3.3.1), a linear fit 
Ao = DqvjT yields Hqw = 10.89 ±0.18, in excellent agreement with Equa¬ 
tion (jlj). The validity of the Gnevyshev-Waldmeier rule is thus also proven 
for individual sunspots. 

Note that from Equation (j^) with the value Co = 32 derived in the 
previous section a mean decay rate of D = 16 would follow. This does not 
contradict to the value of Hgw derived here, as here we are dealing with the 
ratio of the maximal area to the full lifetime of the spot which is typically 
larger than the decay time by about 40 % (Vitinsky, Kopecky, and Kuklin, 
1986). 

The data obviously show a rather large scatter around the linear law, 
and the presence of exclusion zones corresponding to lifetimes around half¬ 
integer rotations makes it even more difficult to discern the correlation, par¬ 
ticularly for samples limited in T —e.g. Vitinsky (1956). However, a mean 
lifetime independent of maximal area, T = constant, is definitely excluded 
by the data ( Q < 1CV 5 , as opposed to Q = 0.18 for the linear best fit). 

Some explanation may be in order as to our choice of Aq as independent 
variable for the fit. Using T as independent variable, as previous authors 
did, would clearly make it possible to avoid the use of a Gnevyshev-Ringnes 
correction, and it would not be sensitive to the short lifetime selection effect 
(Section 3.1): this latter effect leads to the slight upwards deviation of the 
leftmost binned data point in Figure ||, and it is the reason why the data 
were cut down at Ao = 30 in the figure. These advantages of a T-dependent 
fit are however offset by the less uniform coverage of the abscissa in that 
case (the low number of recurrent spots would lend a very small weight 
to those points in the fit) as well as by the possible strong influence of 
the visibility reduction of lifetime (Section 3.1) for some T values (mainly 
for those just below the gap). This influence may explain the suggestion 
by Alexander (1944) and Ringnes (1964a) that groups with a lifetime of 8 
days significantly depart from the linear law. In our material we also find 
somewhat similar deviations near the exclusion zone, but the corresponding 
large error bars indicate that these are not significant. 


6. Conclusion 

In our statistical study of the decay of individual sunspots based on DPR 
data for the years 1977-78 we found that the instantaneous area decay rate 
is related to the spot radius r and the maximal radius ro as 

D = C D r/r 0 C D = 32.0 ± 0.26. (23) 
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This implies that sunspots on the mean follow a parabolic decay law; the 
traditional linear decay law is clearly excluded by the data. A “universal” 
decay law of the form D oc r is also excluded. The validity of the Gnevyshev— 
Waldmeier relationship between the maximal area Aq and lifetime T of a 
spot group 


A) = D G wT T G w — 10 (24) 

is also demonstrated for individual sunspots. No supporting evidence is 
found for a supposed supergranular “quantization” of sunspot areas. 

The robustness of these findings is worth stressing. All the subtleties of 
the data processing described in Section 3 will only lead to slight quantita¬ 
tive modifications in the correlations which are already present in the raw 
data. 

Comparing these results with the model predictions of Table I we find 
that the turbulent erosion model is compatible with all the observational 
data while all the other models are excluded. The turbulent erosion model 
(Petrovay and Moreno-Insertis, 1997) assumes a strict cylindrical geometry 
with no dependence of the field on the coordinate along the tube. This is 
clearly rather far from being a good representation of real sunspots. Yet, 
the fact that the models are able to predict many of the essential features 
of sunspots, including the spontaneous formation of a current sheet around 
the spot and the qualitative and quantitative characteristics of the decay, 
suggests that the model may correctly grasp the essential underlying physics. 

In this paper we have dealt with the mean relationships governing the 
decay of sunspots. In each case however a significant physical scatter is also 
present in the data. While this scatter is apparently “random”, it is already 
known that the decay properties of sunspots may correlate with their other 
physical characteristics. Specifically, a higher decay rate has been found to 
be associated with an irregular shape (Robinson and Boice, 1982), /-type 
polarity (Royal Greenwich Observatory, 1925), more bright structures in 
the umbra (Zwaan, 1968), or higher proper motion (Howard, 1992). The 
problem of random and systematic deviations from the mean laws will be 
treated in the following papers of this series. 
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